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(57) ABSTRACT 

Recently, a smart antenna, i.e., a blind adaptive antenna 
array, has attracted much attention to suppress multiple 
access interference (MAI) and multipath signals to improve 
the capacity of a code division multiple access (CDMA) 
wireless communications system. Most of the smart antenna 
algorithms either need matrix computation or complex cal- 
culations. In this invented direction-of-arrival (DOA) track- 
ing algorithm only 2M comple*x multiplications are required 
per snapshot. Further more the DOA tracking algorithm is 
not sensitive to the mixer phase distortions. Hence the 
system complexity is reduced since a separate phase cali- 
bration is not required for the DOA tracking algorithm. Also 
the equivalent DOA tracking error due to the mixer phase 
distortions is derived. Simulation results show that the DOA 
tracking algorithm in the present invention works effectively 
under additive white Gaussian (AWGN) as well as Rayleigh 
fading environments. Also the system employing the 
invented smart antenna has the same bit error rate (BER) 
performance as those of other existing smart antennas, 
however, has less computation burdens and complexity. 
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DIRECTION OF ARRIVAL ANGLE of the present invention is not sensitive to mixer phase 

TRACKING ALGORITHM FOR SMART distortions. Thus, a separate phase calibration is not required 

ANTENNAS in the system. The proposed algorithm also has lower 

" computational burdens, as compared to conventional 

BACKGROUND OF THE INVENTION 5 algorithms, for calculating weight vectors for smart antenna 

1. Field of the Invention arrays, thereby lowering the complexity and cost of a 
The present invention relates to wireless telecommunica- system. 

tions in general and more particularly to a method and rrtff DF*?Pi?TPTinM hp tup nDAWiwrQ 

system for determining a direction of arrival angle with a BRIEF DESCRIP ™N OF THE DRAWINGS 

smart antenna in a wireless telecommunications system. 10 The aforementioned advantages and features of the 

2. Related Art present invention will be more readily understood with 
A smart antenna, such as a blind adaptive antenna array, reference to the following detailed description and the 

is a method and system used to enhance a desired user accompanying drawings in which: 

signal, to suppress multiple access interfering signals, and in 15 FIG. 1 is a block diagram of a complex pseudo-noise (PN) 

general to improve the performance of a wireless commu- spread receiver for a CDMA communication system con- 

nications system, including CDMA systems such as third figured in accordance with one embodiment of the present 

generation (3G) wireless communication systems, particu- invention. 

larly those specified in the wide band CDMA (W-CDM A) or FIG. 2 is a flow chart of a proposed adaptive DOA 

the CDMA2000 standards. In such systems, smart antenna 20 tracking algorithm in accordance with one embodiment,of 

technology may optionally be employed at a base station. A the present invention. 

smart antenna can suppress interfering signals of different nG . 3is!l ^ pattern £or an M _ 5 elements ^tenna 

uuas ^uireclions ot Arrival; by using spabal diversity. m and DOA=fJ . configured 

in accordance with 

Smart antenna technologies attract much attention these Qne embodiment of the present inve ntion. 

days as they support more users with a high quality of ? r . . , e w ^ , 

service and high data rates. Efficient smart antenna schemes FIG ' 4 15 a *f ™ P a "f ™ for f M=2 elemenlS antenna 

have appeared recently. It is desirable that a smart antenna m ^ and DOA-60 . configured in accordance with 

be effectively easy to implement. one embodiment of the present invention. 

The DOA tracking capability, convergence, and compu- FIG ' 5 is a ^P 0 of the Probabilities of correct direction 

tation complexity are very important issues in the design of 30 detectlon at thc imtial vcrsus thc bit-energy-to-noise 

the smart antennas. Conventional methods for determining P owcr spectrum density E^ in dB for antenna arrays of 2, 

the direction of arrival of an incoming signal in a smart 4 > and 6 elements, configured in accordance with one 

antenna system require matrix computation or complex embodiment of the present invention, 

computation, which are the costly parts of the smart antenna FIG. 6 shows simulation DOA tracking results versus 

tracking algorithms. Specifically, known algorithms attempt 35 iteration with averaging over 100 runs under an additive 

to calculate the weight vectors directly by minimizing or WQ ite Gaussian noise (AWGN) channel, 

maximizing some cost criteria such as SINR (signal to FIG. 7 shows simulation DOA tracking results versus 

interference+noise ratio), BER (bit error rate) or MSE (mean iteration with averaging over 100 runs under a Rayleigh 

square error). These methods require determining a separate fading environment for one embodiment of the present 

weight vector for each antenna element; e.g., if there are six 40 invention. 

antenna elements, six unknown variables must be deter- piG. 8 shows simulation BER results versus EJNq under 

mined. Another drawback associated with known methods is an AWGN environment for a smart antenna configured in 

that practical embodiments of these methods require cali- accordance with one embodiment of the present invention, 

bration of mixer phase distortions. cir n , , „™ l4 „_ T 

r FIG. 9 shows simulation BER results versus EJNq under 

SUMMARY OF THE INVENTION 45 a Rayleigh fading environment for a smart antenna config- 

_ , ured in accordance with one embodiment of the present 

The aforementioned needs are addressed to a great extent invention 
by the present invention, which provides a method and 

system for performing an initial DOA estimate and fine- DETAILED DESCRIPTION OF THE 

tuning the estimate and for tracking the DOA when there is 50 PREFERRED EMBODIMENT 
relative motion between the transmitter and receiver. The 

present invention utilizes the notion that the maximum In tne following detailed description, numerous details, 

antenna gain can be achieved when the weight vector for the SUCD ^ numbers and orientation of antennas in antenna 

antenna array is matched to the array response vector. In arrays, are set forth in order to provide a more thorough 

highly preferred embodiments, rather than attempting to 55 understanding of the present invention. These details should 

calculate the weight vectors directly as is done in conven- not °e understood to limit the invention, 

tional algorithms, the present invention estimates a DOA The present invention is believed to be particularly useful 

angle and calculates the weight vectors as a function of the in CDMA (code division multiple access) systems and thus 

estimated DOA angle. The estimated DOA angle is then will be discussed in that context below. However, the 

adaptively updated to optimize correlation output, e.g., to 60 invention is applicable to other spread spectrum systems and 

achieve the maximum output power or maximum signal-to- is therefore not limited to CDMA systems. The following 

interference-plus-noise-ratio (SINR) for the desired user. description discusses a linear antenna array applied to a base 

This method may be used during both the initial search for station, however, the invention is equally applicable to the 

the DOA, in which a coarse adjustment is made to the mobile side. The present invention is applicable to any smart 

estimated DOA angle, and for tracking the DOA as relative 65 antenna system regardless of the shape of the antenna array, 

movement between the transmitter and the receiver occurs, FIG. 1 is an example receiver structure for a CDMA 

in which a fine adjustment is made. The proposed algorithm communication system with multiple user accesses. A linear 
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antenoa array with M omnidirectional elements is consid- 
ered in this example. The transmitted signal from user k, /2jtdm w t8J 
using the baseband notation, can be written as **(0 = J y B|(r)M0»i«pf->f — — sin^ll* 



where P*', <(>*, f c » a *(0> ^ WO are me transmitted power, j j 2ndm \\ 

phase, carrier frequency, PN chip sequence, and information "W* + 6t " — ^ jj + n " (,) 

data bit sequence for user k, respectively. The phase § k is to 
slowly changing, compared to bit b^t), and time variable for 

<(>* is dropped without loss of generality. Also, the snapshot a nd the complex baseband AWGN -njt) is given by 
rate for weight vector update and the change rate of the DOA 
angle 8 are slower than the bit and chip rate. Therefore, time 

indices for the snapshot, weight vector, and DOA will be 15 ^('^.cW+KjW- ( 9 ) 

also dropped. The PN chip sequence can be written as 

Each antenna element output r m (t) after carrier down con- 
» (2) version (102) in FIG. 1 is multiplied by a complex weight 

MO = 2j a 5 / Vc<'-/ r c) coefficient w m , m«0,l, . . . ,M-1 (103) and then summed to 

20 P r °duce the antenna array output y(t) (104), by assuming 
' that an ideal low pass filter is applied to remove the double 
where a/ is the PN sequence for the W h user and P re (t) is a frequency term. Zero phase distortions at the mixers are 
unit rectangular pulse of chip time interval T c written as assumed for simplicity. Analysis including the mixer phase 

distortions will be presented at the end of the present 
0srsr c (3) 25 invention description. The antenna array output can be 

Pt < U) = \c> othenvisl written as 



M-l (10) 

The information data signal can be written as = itV 2L= ^ w m*M 



30 



and 



where b* is the 1-th data bit for user k and P ri (t) is a unit 
rectangular pulse of bit time interval T b written as 



IN** (") 

XAx Qi x u . . . , x M _A T (12) 

*V r > - \ n otherwise where H and T are the conjugate transpose (i.e., Hermitian) 

40 and transpose, respectively. Assume perfect PN code syn- 
chronization for simplicity. When a weight vector is 
Let k-1 denote the desired user signal index and all radio optimized, the antenna array will yield a beam pattern 
frequency mixers have the same characteristics at the towards the desired user signal and suppress interfering 
receiver. Then, the complex signal received at the m** 45 signals of different arrival angles. 

antenna element, m-0,1, . . . , M-l, can be written as j^e weight vector of the linear antenna array is normal- 

ized to have a unit norm, i.e., ^ W W-L When the antenna 
r m (r) = V2PTai(0fci (0ffiexpf/2*/ e r + tp x - ?^an 9i )] 4. ^ array weight vector is perfect for the desired user signal, the 

weight vector can be written as 



JlP^a k {i - r k )b k (/ - r, )a k 
exp^2ff/ c f + <p k + 6 k - ^HzmQ k jj + nJO 



50 



1 r ■ a .IxdtM-l) . . i r (13) 



The antenna array response vector a is the same as the 

where m-0 is the reference antenna element (101) in FIG. 1, perfect weight vector in equation (13). Because the antenna 

d is the distance between the two adjacent antenna elements, elements are in a linear array with space d, the antenna array 

X is the wavelength of the carrier frequency, a*, and x k output y(t) would be 
are the amplitude, phase, and delay for the k A user, 

respectively, and the bandpass noise n^t) at the m 4 * element 60 < ( i 4) 

can be written as >(0 = %! i x.U) = J ^p- a t (i)fr, (0<*i + #0 + yit) 

njt)-n„^)cce Txfj-n^in 2xfj. (7) 

Since coherent demodulation is assumed for the desired user 

1, 6j is set to zero. 65 where |(t) and v(t) are the interference signals due to the 

The down converted signal at the m** antenna element is other users and the noise component, respectively, and can 

denoted by x m (t) and can be written as be written as 
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m-i r k i — 



(15) 



It is observed from equation (14) that when the weight 
vector is perfectly matched with the array response vector in 
the right hand side of equation (13), the receiver with a linear 
antenna array will have an amplitude gain of M 1/2 , i.e., the 
received signal power with M antenna elements will be M 
times stronger than that with a single antenna receiver. For 
other kinds of antenna arrays, e.g., circular or square arrays, 
etc., the same maximum antenna gain can be achieved when 
the weight vector is perfectly matched with the array 
response vector. The present invention can be applied to the 
other antenna arrays without loss of generality. The smart 
antenna output y(t) is fed to the successive post processor 
(105), e.g., demodulator, to retrieve the user data bit infor- 
mation. 

We observe that y(t) is a Gaussian random process and the 
mean of y(t) is 

EiyiO) = Z + > W')l«p(j ^^sin0,]J (17) 

= ^= X t^W(O) +> ^^(DJlex^ ^^sin^, j = 0, (18) 



and the variance of y(t) is given as 

l ff w ~ l 
£l7<0/('» = -£| £ow 



/ 2ndm „ \\ (l9 > 
(0 +yn m ^(r))cxply — ^— sinft || x 



■2* 



(20) 



weight vector W-[w 0 , w 1 



VM 1 J 



(21) 



(16) 



Then the antenna array output signal J(t) (104) can be 
written as 



10 



HO = £ *m*mU) = 5(0 + MM0 + N(f ) 
M-l , 

Z\ P, f tTndm x Txdm \\ 

M-I K I 

MAw= ZZv^ a,< ' ,^,i( ' , 

20 

( f t Txdm . ft 2ndm . „ \\ 
«*«PUlWk +4 + — — %m9i — siiA II 

" (0 = 7=|^<M^H 



(22) 



(23) 



(24) 



(25) 



where S(t), MA(t), and N(t) represent the signal, multiple 
3Q access interference, and thermal noise component, respec- 
tively. The output J(t) is despread with the desired user's PN 
sequence a^t). The signal component of the despread output 
sampled for the X th bit can be written as 



35 



-r 

JiT b 



40 



Pi 1 -<?~ 



(27) 



45 



where 



where a^-var^^OJ-varfn^Q^-No^ is the variance of 
the AWGN thermal noise. 

Most of the existing smart antenna algorithms update the 



Jir, 



a x {i)a v {t)dt=T b . 



(28) 



50 



^m-iY so toat tne weight 



vector minimizes or maximizes a cost or criteria function, 
e.g., minimizes a mean-square-error (MSE), minimizes a bit 
error rate (BER), or maximizes SINR. There are M unknown 
variables in the weight vector \¥. Instead of trying to find the 
M unknown variables, the present invention updates the 
DOA 8 a for user 1 and simply uses the estimate 6j in 
equation (13) to calculate a weight vector. Only one variable 
8j is updated instead of M variables. The present invention 
presents a method to find the optimum estimate 8j adap- 
tively to maximize correlation output with the tentative 
decision feedback bit. Two different step sizes A8 will be 
employed for the updating in the preferred embodiment. For 
the initial search, a large step size will be employed, and for 
the tracking a fine step size will be employed. Thus, the 
weight vector (107) using an estimated DOA 9j for user 1 
can be written as 



When the adaptive estimation 8j approaches 8 1( the 
desired signal component S(1T A ) for the X th bit can be written 
as 



S(lT b ) = 



lim 



(29) 



60 



MP, 



65 



The interference and thermal noise component of the 
despread output sampled for the 1** bit can be written as 
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8 



MA(rr b ) 



Jir b 

K 

Z 



MA(t)a l (t)dt = 



and 



= — f 
\[m~ Jrr. 



N(tMt)dt = 



(31) 



where 



Hit 



a k {t~ 



n) ai (t)b t (T)dt. 



The real part of the despread output at the sampling time 
instance (l+l)T fc is preferably used instead of the magnitude. 
Let Y1(1T 6 ) denote the real part of the weighted and 
despread correlation output. Then, 



K(/r 6 )=/?«{5C/r t> )+Ai4(/r 6 )+JV(rr A )}. 
The MA(lT fc ) and N(lT fr ) have zero means, 



Let Yl denote the expectation of the real part of the 
despread output multiplied by the decision feedback bit b, 3 , 
Then. Yl can be written as 



Yi(Q u Qd-EiYlitTjb, 1 } (36) 

-£{/&{s(/r A )5 / 1 +AM(/r i> )6/^(/r t )fc / 1 }} (37) 

-ReiSitTJbt 1 } (38) 

n(e l j;) = E{Yim)b]} m 

= E {Kc{s{tT b )b l t + MA{iT b jb) + N{lT b )h\)) (37) 

= Re{S(/n$} (38) 

L (40) 



If a pilot channel (106) is available in a CDMA system, 
(30) the pilot channel signal may be used as the desired signal 
may be used to track the DOA of the desired user. The pilot 
channel signal amplitude may flunctuate under fading envi- 
5 ronment. This fading environment effect will be discussed 
below. If a pilot symbol pattern is available, the pilot pattern 
b/, which is known to both the transmitter and the receiver, 
may be used; or the temporary decision feedback informa- 
tion bit b/ may be used for b/. The decision feedback 
to information bit b, 1 is the postprocessor (105) output, which 
is the binary demodulation output in the present invention. 

FIG. 2 shows an example flow chart of the proposed 
adaptive DOA tracking algorithm which is used at the smart 
antenna weight vector update (107) in FIG. 1. Here the 
15 assumption is that the angle is changing slowly, compared to 
the bit rate. This is a reasonable assumption. For example, 
suppose a mobile station is 1 km away from a base station, 
its moving speed is 100 km/hour, the data rate is 9600 bits 
(32) p er second (bps), and the DOA changes 1.5 degrees. Then 
20 (15*Ji-103-3600-9600)/(180100103)=9047 bits will have 
been transmitted during the 1.5 degree change of DOA. For 
higher data rate, more data bits can be transmitted. Suppose 
the weight vector is updated every 20 bits. Then, 450 
updates can be made for 1.5 degree DOA change. We also 
25 assume zero angle spread for each multipath for simplicity. 
If the DOA of the desired signal is known, then the weight 
vector can be automatically obtained by using, for example, 
equation (13). However, the DOA of the desired signal is 
usually changing due to the mobile movement. Hence an 
30 adaptive DOA tracking algorithm is desirable in practice. 
Assume a power control is used so that the received power 
from any user is approximately equal to the others, which is 
a reasonable assumption in practice. The antenna array 
output will have a peak output when the beam is pointing to 
35 the desired user and the PN spread processing gain is larger 
than the number of users. 
The estimation 9 a is updated as 



(33) 

(34) 
(35) 



(41). 




if b/V-L Yl(8 1 ,6 / ) is the function of (sin Oj-sin B,). FIGS. 
3 and 4 show the shape of Yl except the amplitude 
duTerence, for a beam pattern of M-6 and M-2 elements, 
respecively, and DOA-0° and DOA-60 0 , respectively. The 
solid line is the corresponding magnitude of EjYIQTJb/}, 
and the dashed line is the real part of EjYlOT^b, 1 }. The 
real part has an smaller beamwidth than that of the magni- 
tude which means that the real part has better resolution on 
the DOA. After an intitial estimation of the DOA from the 
desired user is found, the DOA tracking algorithm such as in 
FIG. 2 is activated to track the desired user. As a result, the 
beam will focus on the desired user to keep the maximum 
correlation output using the antenna output. 



40 where i denotes the update index. It is not nessessary to track 
the DOA at chip or bit rate since^ the DOA changes slowly, 
compared to the bit rate. When 8^-9,, the weight vector is 
steered pointing to the desired user. 

Preferred embodiments of the initial DOA estimation 

45 process (also referred to as the DOA acquisition process) 
and the DOA tracking process will now be discussed in 
further detail. As discussed above, the initial DOA estima- 
tion process is necessary when the base station does not 
know the DOA angle for a desired user. The goal of the DOA 

50 estimation process is to estimate the DOA such that the 
estimate falls within the main lobe of the transmitted signal. 
Generally, the process comprises the following steps: form 
an initial estimate of the DOA angle, calculate a weight 
vector using this initial estimate, obtain the correlator output 

55 over some period using the weight vector, and repeat the 
previous steps until the estimated DOA is acceptable (i.e., 
within the main lobe of the transmitted signal). 

In one highly preferred embodiment, the estimated DOA 
angle is initialized to an angle corresponding to the sector 

60 boundary and then stepped (incremented or decremented, 
depending upon the initial estimate) by a coarse step AO until 
the entire sector has been traversed. It should be noted that 
what is meant by "an angle corresponding to a sector 
boundary" includes both an angle equal to the sector bound- 

65 ary as well as an angle approximately equal to the sector 
boundary, such as an angle that is within some fraction (e.g., 
one half) of a coarse step A8 of the sector boundary. The 



05/05/2004, EAST Version: 1.4.1 



US 6,483,459 Bl 



10 



coarse step increment AO is preferably set to one half of the 
beamwidth of the main lobe, as shown in the fifth column of 
Table 1 below. This is done to ensure that at least one of the 
initial estimates is within the main lobe such that the 
correlator output will exceed the correlator output for other 
angles corresponding to side lobes. The initial estimate is 
used to obtain a trial set of weight vectors. The correlator 
output is then obtained using the trial weight vectors over 
one bit (symbol) period. The estimated DOA angle is then 
updated, preferably every bit time interval, with a coarse 
step AG. This process continues until the entire sector has 
been traversed. The angle with the largest correlator output 
is then chosen as the acquired DOA and the tracking 
algorithm begins. The number of steps to search 180 degrees 
(it is not necessary to search through 180 degrees where a 
sector size is smaller, such as 120 degrees) is shown in Table 
1 for different numbers of antenna elements. At most only 
about 25 iterations (25 bit time intervals) are required to 
bring the beam pattern pointing to the desired user from any 
worst initial angle. This method is especially effective in 
systems in which power control is not used. 

Alternatively, rather than comparing correlator outputs for 
several DOA angles, comparing them to each other to 
choose the corresponding DOA angle, correlator outputs for 



10 



elements is used. In other embodiments, the fine step incre- 
ment A8 is based on the angle error. In one highly preferred 
embodiment, the fine step increment is set to 10% of the 
angle error from Table 1 corresponding to the number of 
elements in the array. 

Reference is made to FIG. 2. The estimated DOA ange is 
initialized for two trial angles at Step 200. The following is 
a preferred updating (tracking) procedure for 8 2 at the i* 
snapshot: 

Step 201: Generate two sets of weight vectors^ Wl and W2 
as shown in FIG. 2 with the estimated Q v Here the 
subscripts 1 and 2 represent two trial angles; early angle 
Bj-Ae/2 and and late angle Oi+AO/2, respectively. 



15 



20 



. M-l 



w2„ = 



i*TJ, m = 0, 1 M-l. 



(42) 



(43) 



Step 202: Measure two despread outputs Yl(lT fc ) and 
Y2(\T b ) in equation (33) at each sampling time instance 
(1+1)T 6 by using ^1 and W2 respectively, and X in 
equation (12) and PN sequence. 



estimated angles could be compared to a threshold. In such g 203 . Take me of YlOTJb, 1 and Y20TJ6, 1 in 



a scheme, the first estimated DOA angle with a correlator 
output exceeding the threshold is selected. For example, a 
first estimate DOA is made (this could be done using a 
random process, or the first estimate could correspond to a 
sector boundary as discussed above) and the corresponding 
correlator compared to the threshold. If the threshold is 
exceeded, that angle forms the initial estimate and the 
tracking process begins. Otherwise, the estimated angle is 
updated again (either randomly or by the coarse step size) 
and the process continues until the threshold is crossed. 

Table 1 shows example antenna array beam pattern 
parameters, for a 180° search range, where the number of 
antenna elements varies from 2 to 10. 

TABLE 1 

Antenna array beam pattern parameters 



equation (36) over L bit time intervals, 



30 



Y2= [ l f j Y2(lT b ji>l 



(44) 



(45) 



35 Step 4: Update the estimated angle Q x as the following 



d,(f+i) = 



40 



0i(O + A0. if YI<Y2 
3i(/)-A0, iff/ > 72 
©i (/), otherwise, 



(46) 



Number 

of 
antenna 
elements 



Amplitude 
Angle Relative 
error to the 
(in maximum 
degrees) gain 



Number of 
steps in 
the 180° 
search range 
during the initial 
search period 



A6 (in degree) 
^ main 
lobe beamwidth x Vi 
during the initial 
search period 



2 


11.0 


91.3% 


3 


60 


3 


6.0 


91.3% 


6 


31.8 


4 


4.4 


90.2% 


9 


21.6 


5 


3.3 


90.5% 


11 


16.4 


6 


2,7 


90.3% 


14 


13.2 


7 


2.3 


90% 


17 


11 


8 


1.9 


90.8% 


19 


9.6 


9 


1.7 


90.5% 


22 


8.4 


10 


1.5 


90.7% 


25 


7.4 



45 



50 



55 



When an acceptable estimated DOA is acquired, the DOA 
tracking algorithm is activated with a fine step increment A8. 
The estimated DOA during the initial acquisition period 
used as the initial value for the fine tracking period. The 60 
choice of fine step increment A0 depends on the number of 
antenna elements. Generally, the larger the number of 
antenna array elements is, the smaller the A8 will be 
desirable. It is a function of the number of array elements M. 
In some embodiments, the fine step increment A9 is set to 2° 65 
because the bit error rate is affected insignificantly by the 
beam direction error of 2° when an antenna array of M-2 



In step (203) the time average is used to approximate the 
ensemble average in equation (36). At step 205, the weight 
W is updated using the new estimated angle. 

The above-described tracking process (steps 201 to 204) 
ideally continues for as long as the incoming signal is 
received. However, under some circumstances, tracking 
may be lost, such as when the DOA may suddenly change 
due to a change in the environment, or the initial DOA 
estimate may correspond to a sidelobe rather than the main 
lobe. It is desirable to detect such events and to re-acquire 
the DOA by triggering the acquisition process. One indica- 
tion that could be used to trigger the DOA acquisition 
process is that the received correlation output after despread 
is less than a threshold, e.g., 50% of the peak gain of the 
antenna array. Another scheme is to trigger the acquisition 
process when the bit error rate or the frame error rate (or 
some other signal quality measurement) crosses a threshold. 
A third scheme is to keep the acquisition process running in 
parallel and reinitialize the tracking process if the correlator 
output for any estimated angle exceeds the correlator output 
for the current angle determined by the tracking process. 

In practical channel environments, the desired user signal 
may come from different DOAs with time delays. If the 
desired signals are from the same DOA with different time 
delays, then a RAKE receiver structure can be employed to 
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collect the resolvable signals. If the desired signals are from 
different DOAs, several duplicate weight vector generators 
(107 in FIG. 1) can be employed to collect the desired user 
signals received from different resovable DOAs. 
Fading 

The proposed DOA tracking algorithm discussed above is 
still valid under Rayleigh fading environment. The reason is 
that the fading amplitude of the desired user only affects the 
amplitude of the despread output at each sampling instance. 
For a given DOA, if Yl(iT fa )<Y2(lT fc ), then a u Yl(lT 6 )<a^ 
/Y2(lT fc ), since a iy >0. Even with a u changing during L bit 
time interval, the following equality still holds: 



12 



Let i|> denote the equivalent DOA tracking error due to the 
mixer phase distortions, and 6, and 6, be the true and the 
estimated DOA, respectively. Then, 9 a =9 j+tJj . The degraded 
output Yl due to the equivalent DOA tracking error can be 
written as 



F/ = 



(54) 



10 



1-0 /-o 



(47) 



15 



sun) 



b 



S{t)aiU)dt 



2M 



m-0 



and 



Yl = £{Re[5(/7;)5;]} 

m=0 



The equivalent DOA error due to the mixer phase distor- 
tions can be determined by the following equation: 



Hence the DOA tracking procedure remains unchanged. 
No Phase Calibrations 

In a practical system, the characteristics of the mixers are 
different and they introduce different phase distortions at the 
mixers in general. Let ip OT be the phase distortion due to the 
m f * mixer. Without loss of generality, uj is assumed to be 
uniformly distributed between [-Q\ +Q ] and {ty m } are 
independent. Including the mixer phase distortions, the 
signal component of the antenna array output after PN 
despreading can be written as 



20 



C*tf|+*>-«fn0i>| Msinn 



(55) 



25 



by using eqution (26). 

When the beam pattern is pointing to the DOA of the 
desired user, the desired signal output of the antenna array 
can be written as 



(51) 



sinfl 

i2M , |£|ail «_ 



The degradation of the gain due to the mixer phase 
distortions is about sin Q/Q. For example, when Q»60°, the 
degradation is only 10 log(0.827)— 0.825 dB. The angle 
error due to the mixer phase distortions depends on the 
number of the antenna array elements. 



where uj is the DOA error due to the mixer phase distortions. 
As a special case, let S^O 0 . Then, the equivalent DOA 
tracking errors due to the mixer phase distortions are listed 
in Table 2 for M=2, 3, 5, 7 antenna elements. 

TABLE 2 



(48) 30 


Equivalent DOA tracking error due to the mixer phase distortions with 




mixer phase difference of ranee Q = 60° 


(49) 


M 






2 


15.85° 


35 


3 


8.6° 




5 


4.5 




7 


3.07° 



40 



(50) 



45 



50 



(52) 



55 



(53) 



The tracking algorithm will work if the initial DOA guess 
is within the main lobe. The main lobe width of the beam 
pattern will affect the coarse acquisition time. The wider the 
main lobe width, the shorter the acquisition time. When the 
DOA is larger than 60° or less than -60°, the initial DOA 
search might take a wrong direction due to the poor reso- 
lution of the antenna array beam patern using the DOA. If 
the true DOA of the desired user is about 60°, a correlation 
output with a DOA guess in [-60 o ,-70 o ] might be higher 
than that with a guess of around 60° due to noise and MAI. 
This can be observed from FIG. 4. Fortunately a typical base 
station employs three sectors. Then, the DOAs for each 
sector can be reasonably restricted to within [-60°, 60°]. 

FIG. 5 shows the probabilities of correct direction detec- 
tion at the initial search versus the bit-energy-to- noise power 
spectral density ratio E^/Hq in dB for antenna arrays of 2, 4, 



60 spectral density in watts per Hertz. For simulation, binary 
phase shift keying (BPSK) modulation is employed. An 
antenna array of M«2 elements is used. A total of nine users 
are assumed by way of example in the system. The data rate 
is 9600 bps. The carrier frequency f c is 1.9 Giga Hertz. A 

65 Gold sequence of length 31 is used as the PN sequence. 
Under Rayleigh fading environment, the mobile speed V is 
set to 50 km/hour. 
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TABLE 3 



List of the number of computations per 
snapshot for different smart antennas 



Adaptive algorithms 



Number of complex 
computations/snapshot 

Multiplication Addition 



2M 
2M 
2M 3 
2M 



+ 2M - 1 

- 1 



10 



15 



20 



25 



30 



Least-Mean-Square (LMS) 2M + 3 

Normalized LMS 3M + 1 

Recursive Least-Square (RLS)* 2M 2 + 3M 

Optimum adaptive convergence 4M + 3 
parameter 

Maximum output power criteria without 3M 4- 1 3M - 1 

Lagrange multiplier 

Present invention** 2M 2M - 2 

D. Shim and S Choi 5M + 1 3M - 1 

blind adaptive algorithm 



•Direct RLS requires about 2M 2 + 5M operations per snapshot 
*"The present invention does not update the DOA every chip or bit. 

Table 3 compares the list of the number of computations 
per snapshot for different smart antenna algorithms where M 
is the number of antenna array elements. Table 3 shows that 
the present invention has a smaller number of computations' 
per snapshot than most of the existing smart antenna algo- 
rithms. 

By updating the DOA every L bits where L=20 bits were 
used throughout the invention, FIGS. 6 and 7 simulation test 
results for show the transient as well as steady state of the 
DOA tracking behavior under AWGN and Rayleigh fading 
environments, respectively, when M=2, Ej/N 0 =6 dB, true 
DOA=30° and initial DOA estimation was 50°. In FIG. 7, 
mobile speed V=50 km/hour, and carrier frequency f c =1.9 
Giga Hertz. We observe that for both AWGN and fading 
environments, the proposed algorithm takes only about 10 35 
and 18 iterations when the mixer phase distortions are absent 
and present, respectively. 

FIGS. 8 and 9 show simulation BER results versus Ej/N 0 
under AWGN and Rayleigh fading environments, 
respectively, by employing the present DOA tracking 
algorithm, maximum output criteria without Lagrange mul- 
tiplier (3), and Shim and Choi's blind adaptive algorithm 
(4). The simulation results indicate that the smart antennas 
configured in accordance with one embodiment of the 
present invention have the same performance, compared to 
the other smart antenna algorithms. In these simulations, 
BPSK modulation is applied. M=2 antenna elements array is 
used. A total nine users are assumed in the system. The data 
rate is 9600 bps. The carrier frequency i c is 1.9 Giga Hertz. 
Walsh sequences of length 32 are used as the PN sequences. 
Under Rayleigh fading environment, the mobile speed V is 
set to 50 km/hour. 

What is claimed is: 

1. A method for acquiring a direction of arrival angle for 
an incoming spread spectrum signal comprising the steps of: 
initializing an estimated direction of arrival angle of an 

incoming spread spectrum signal; 
calculating a weight vector for a smart antenna array 

based on the estimated direction of arrival angle, the 

weight vector being matched to an array response 

vector for the antenna array; 
weighting signals received by the smart antenna array 

with the weight vector to form a weighted output 

without performing phase calibration; 
despreading the weighted output to form a despread 

output; and 



45 



50 
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selecting the estimated direction of arrival angle as a 
direction of arrival angle based upon the despread 
output. 

2. The method of claim 1, wherein the selecting step is 
performed by comparing the despread output to a threshold. 

3. The method of claim 2, further comprising the steps of: 
updating the estimated direction of arrival angle if the 

estimated direction of arrival angle is not selected in the 
selecting step; and 
repeating the updating, calculating, weighting, and 
despreading steps if the estimated direction of arrival 
angle is not selected in the selecting step'. 

4. The method of claim 3, wherein the estimated direction 
of arrival angle is changed by a step size. 

5. The method of claim 4, wherein the step size is 
approximately one half of a beamwidth of the incoming 
spread spectrum signal. 

6. The method of claim 3, wherein the estimated direction 
of arrival angle is modified by randomly selecting a new 
direction of arrival angle within a sector corresponding to 
the smart antenna array. 

7. The method of claim 1, wherein the estimated direction 
of arrival angle is initialized to correspond to a boundary of 
a sector in the initializing step. 

8. The method of claim 7, further comprising the steps of 
updating the estimated direction of arrival angle by a step 
angle and repeating the updating, calculating, weighting, 
and despreading, steps until the estimated direction of 
arrival angle has traversed the sector. 

9. The method of claim 8, wherein the step angle is 
approximately one half of a beamwidth of the incoming 
spread spectrum signal. 

10. The method of claim 8, wherein the selecting step is 
performed by selecting the estimated direction of arrival 
angle corresponding to the largest despread output from 
among all estimated direction of arrival angles in the sector. 

11. The method of claim 1, wherein the array is linear. 

12. The method of claim 11, wherein each element of the 
weight vector is calculated according to an equation of the 
form 



40 



where W v is a weight for Mth element of the array, C is a 
constant, M is an antenna array element number, X is a 
wavelength of the incoming spread spectrum signal, d is an 
antenna array element spacing, and 6^ is the estimated 
direction of arrival angle. 

13. The method of claim 11, wherein the weight vector is 
calculated according to an equation of the form 



55 



wherein W is the weight vector, M is an antenna array 
element number, X is a wavelength of the incoming spread 
spectrum signal, d is an antenna array element spacing, 0 A is 
the estimated direction of arrival angle, and T is the trans- 
pose. 

14. The method of claim 1, wherein an imaginary portion 
of the despread output is not used in the selecting step. 

15. The method of claim 1, wherein both a real and an 
imaginary portion of the despread output is used in the 
selecting step. 

16. The method of claim 1, wherein the incoming spread 
spectrum signal is a code division multiple access signal. 
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17. The method of claim 1, wherein the incoming spread 28. The method of claim 27, wherein each element of the 
spectrum signal is a pilot channel signal. first weight vector is calculated according to an equation of 

18. The method of claim 17, wherein the despread output the form 
is multiplied by a pilot symbol pattern. 

19. The method of claim 1, wherein the incoming spread 5 _ tW 
spectrum signal is multiplied by a tentative decision feed- ^ M ~ Ce 

back bit. 

20. The method of claim 1, wherein the signals weighted wn cre is a weight for Mth element of the array, C is a 
in the weighting step are collected over a bit period. constant, M is an antenna array element number, X is a 

21. A method for tracking a direction of arrival angle to wavelength of the incoming spread spectrum signal, d is an 
comprising the steps of: antenna array element spacing, and 6, is the first trial angle. 

initializing an estimated direction of arrival angle of an 29. The method of claim 27, wherein the first weight 

incoming spread spectrum signal; vector is calculated according to an equation of the form 

forming a first trial angle, the first trial angle being spaced 

apart from the estimated direction of arrival angle by a 15 w _ 1 r jlj^, j *"**- 1 ' ^, i r 

first angular difference in a first direction; ~~ ~ V£T ^ ' * * * 

calculating a first weight vector for a smart antenna array 

based on the first trial angle, the first weight vector whereia w ^ ^ firs , we ^, M is an antenna array 

being matched to an array response vector for the ek - mber> x k a wa Velength of the incoming spread 

antenna arrav or 

J ' spectrum signal, d is an antenna array element spacing, 8j is 

weighting signals received by the smart antenna array the first trial angle, and T is the transpose, 

with the first weight vector to form a first weighted 30. The method of claim 21, wherein an imaginary portion 

output without performing phase calibration; 0 f the first despread output is not used in the selecting step. 

desp reading the first weighted output to form a first 25 31. The method of claim 21, wherein both a real and an 

despread output; and imaginary portion of the first despread output is used in the 

selecting the first trial angle as a new direction of arrival selecting step. 

angle based upon the first despread output. 32 - The method of claim 21, wherein the incoming spread 

22. The method of claim 21, further comprising the steps spectrum signal is a code division multiple access signal, 
of: 30 33. The method of claim 21, wherein the incoming spread 

forming a second trial angle, the second trial angle being s ^?™ signal » a /* lot c !j aDnel si ^ al * 
spaced apart from the estimated direction of arrival . 34 m meu^od of claim 33 wherem the despread output 
angle by a second angular difference in a second 15 ™ l !£ hGd a V 1 } 0 } s ^° l P attera . 
direction, the second direction being different from the 35 ' ^ method of claim 21 > wherein the incoming spread 
first direction* 35 sp^ 1 ™ 03 signal is multiplied by a tentative decision feed- 
calculating a second weight vector for a smart antenna ^jj* .i_ j * i • L • . r . • 
„ , ° , 4 , j. ■ i , j ■ Li 36. The method of claim 24, wherein a loss of tracking is 
array based on the second trial ancle, the second weight j * . j i_ * L ./ ^ L ■ - 

/ , . 4 i i , , % detected by comparing a bit error rate for the incoming 

vector being matched to an array response vector for , . . f. *- i_ u 

the antenna arra * spread spectrum signal to a threshold, 

e an enna array, ^ ^7 -j^g mem od of claim 24, wherein a loss of tracking is 

weighting signals received by the smart antenna array detected by ^m^ting a frame error rate for the incoming 

with the second weight vector to form a second spread spectrum signal t0 a threshold, 

weighted output; and 3g ^ method of daim 24 wherein a loss of trackin g ; s 

despreading the second weighted output to form a second detected by comparing a signal strength to a threshold. 

despread output; 45 39. A method for acquiring a direction of arrival angle for 
wherein the selecting step is performed by selecting the an incoming spread spectrum signal comprising the steps of: 
first trial angle as a new estimated direction of arrival initializing an estimated direction of arrival angle of an 
angle if the first despread output is larger than the incoming code division multiple access spread spec- 
second despread output, selecting the second trial angle signa i ta corre spond to a boundary of a sector 
as a new estimated direction of arrival angle if the so calculatin a wei ^ t vector for a smart antenna 
second despread output is larger than the first despread based Q * ^ directioQ of ^ ^ 
outpu , and otherwise leaving the estimated direction of wei ^ t vectof bei matched aQ ^ 

arnvalar^le unchanged vector for the antenna array; 

23. The method of claim 21, further comprising the steps ...... . 

of repeating the forming, calculating, weighting, despread- 55 wei ghting signals received by the smart antenna array 

ing and selecting steps until a loss of tracking is detected. over a s y mbo1 V* no6 Wlth the wei S ht vector t0 form a 

24. The method of claim 22, further comprising the steps weighted output without performing phase calibration; 
of repeating the forming, calculating, weighting, despread- decorrelating the weighted output to form a decorrelated 
ing and selecting steps until a loss of tracking is detected. output; 

25. The method of claim 21, wherein the selecting step is 60 multiplying the decorrelated output by a symbol bit corre- 
performed by comparing the first despread output to a sponding to a symbol represented by the incoming code 
threshold. division multiple access spread spectrum signal during the 

26. The method of claim 21, wherein the first angular symbol period to form a despread output; 

difference is approximately equal to ten percent of an angle changing the estimated direction of arrival angle by a step 

error of the smart antenna array. 65 size approximately equal to one half of a beamwidth of 

27. The method of claim 21, wherein the smart antenna the incoming code division multiple access spread 
array is linear. spectrum signal and repeating the calculating, 
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weighting, decorrelating, and multiplying steps for 
each estimated direction of arrival angle until the sector 
has been traversed; 
selecting the estimated direction of arrival angle corre- 
sponding to the largest despread output as the direction 5 
of arrival angle. 

40. The method of claim 39, wherein the array is linear 
and the weight vector is calculated according to an equation 
of the form 

10 

wherein ¥L is the weight vector, M is an antenna array 
element number, X is a wavelength of the incoming spread 15 
spectrum signal, d is an antenaa array element spacing, 6 3 is 
the estimated direction of arrival angle, and T is the trans- 
pose. 

41. The method of claim 40, wherein only a real portion 

of the decorrelated output is multiplied by the symbol bit to 20 
form the despread output. 

42. The method of claim 41, wherein the incoming spread 
spectrum signal is a pilot channel signal and the symbol bit 
corresponds to a pilot symbol pattern. 

43. The method of claim 41, wherein the symbol bit 25 
corresponds to a tentative decision feedback bit. 

44. A method for tracking a direction of arrival angle 
comprising the steps of: 

initializing an estimated direction of arrival angle of an 
incoming spread spectrum signal to a value within a 30 
main lobe of the incoming spread spectrum signal; 

forming a first trial angle, the first trial angle being spaced 
apart from the estimated direction of arrival angle by a 
step angle in a first direction; 

forming a second trial angle, the second trial angle being 
spaced apart from the estimated direction of arrival 
angle by the step size in a second direction; 

calculating first and second weight vectors for a smart 
antenna array based on the first and second trial angles, ^ 
respectively, the weight vectors being matched to an, 
array response vector for the antenna array; 

weighting signals received by the smart antenna array 
during a symbol period with the first and second weight 
vectors to form a first and second weighted outputs 45 
without performing phase calibration; 

decorrelating the first and second weighted outputs to 
form first and second decorrelated outputs; 

multiplying a real portion of the first and second decor- 
related outputs by a symbol bit corresponding to a 50 
symbol represented by the incoming code division 
multiple access spread spectrum signal during the sym- 
bol period to form first and second despread outputs; 

selecting the first trial angle as a new estimated direction 
of arrival angle if the first despread output is larger than 55 
the second despread output, selecting the second trial 
angle as a new estimated direction of arrival angle if the 
second despread output is larger than the first despread 
output, and otherwise leaving the estimated direction of 
arrival angle unchanged; and 
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repeating the calculating, weighting, decorrelating, 
multiplying, and selected steps until a loss of tracking 
is detected. 

45. A method for tracking a direction of arrival angle for 
a signal comprising the steps of: 

estimating an estimated direction of arrival angle; 
calculating a weight vector for an antenna array compris- 
' ing a plurality of antenna elements by matching the 

weight vector to an antenna array response vector based 

on the estimated direction of arrival angle; 
obtaining a signal weighted by the weight vector from the 

antenna array without performing phase calibration; 
updating the estimated direction of arrival angle based on 

a parameter of the signal. 

46. An apparatus for acquiring a direction of arrival angle 
for an incoming spread spectrum signal, the apparatus 
comprising: 

a processor connectable to a smart antenna array, the 
processor being configured to perform the steps of 
initializing an estimated direction of arrival angle of an 

incoming spread spectrum signal; 
calculating a weight vector for the smart antenna array 
based on the estimated direction of arrival angle, the 
weight vector being matched to an array response 
vector for the antenna array; weighting signals 
received by the smart antenna array with the weight 
vector to form a weighted output without performing 
phase calibration; 
despreading the weighted output to form a despread 
output; and 

selecting the estimated direction of arrival angle as a 
direction of arrival angle based upon the despread 
output. 

47. The apparatus of claim 46, further comprising a smart 
antenna array connected to the processor. 

48. An apparatus for tracking a direction of arrival angle, 
the apparatus comprising: 

a processor connectable to a smart antenna array, the 
processor being configured to perform the steps of 
initializing an estimated direction of arrival angle of an 

incoming spread spectrum signal; 
forming a first trial angle, the first trial angle being 
spaced apart from the estimated direction of arrival 
angle by a first angular difference in a first direction; 
calculating a first weight vector for a smart antenna 
array based on the first trial angle, the first weight 
vector being matched to an array response vector for 
the antenna array; 
weighting signals received by the smart antenna array 
with the first weight vector to form a first weighted 
output without performing phase calibration; 
despreading the first weighted output to form a first 
despread output; and selecting the first trial angle as 
a new direction of arrival angle based upon the first 
despread output. 
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